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Abstract 

We have studied the behavior of the 5-matrix poles near threshold for quantum waveguides 
coupled to a cavity with a defect. We emphasize the occurrence of both dominant and shadow 
poles on the various sheets of the energy Riemann surface, and show that the changes of the total 
conductivity near threshold as the cavity's width changes can be explained in terms of dominant 
to shadow pole transitions. 

PACS numbers: 73.63.Nm Quantum wires, 73. 23. Ad Ballistic transport, 72. 10. Fk Scattering by point 
defects, dislocations, surfaces, and other imperfections (including Kondo effect) 
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In multichannel scattering, the poles of the S-matrix lie in general on different sheets of 
the Riemann energy surface. Poles on the negative real axis of the physical sheet correspond 
to bound states of the system, while poles on the unphysical sheets close to the physical 
energy axis are associated to resonant states. In addition to the resonance, dominant poles, 
however, there are also poles on unphysical sheets, far away from the physical region, which 
do not have in general observable effects, and have been referred to as shadow poles |l[ . The 
possible role of shadow poles has been discussed over the years in the context of particle 
2], nuclear [2I and atomic [4] physics, as well as in laser-induced multiphoton processes 
5|, |6j . Their study is of particular relevance when the scattering process depends upon 
some tunable parameter; as this parameter is changed, the S'-matrix poles move on the 
various sheets of the energy Riemann surface, and may pass a scattering threshold. In so 
doing, some shadow pole may approach the physical region, thereby becoming dominant, 
and producing observable effects, whereas a previously dominant pole may retire to a less 
exposed position Q. 

In this paper we would like to point out another situation, where shadow and dominant 
poles may exchange their roles, and give rise to non-trivial observable effects near threshold. 
It stems from recent developments in nanotechnology, which allow one to obtain a strictly 
two-dimensional electron gas subject to confined geometries js[ To be definite, we shall 
consider the device of Fig. [I], where a resonant cavity or stub having width c and length l s 
is coupled to a uniform guide of indefinite length and width b. The stub contains a region, 
depicted by the shaded area in Fig. [IJ with a defect described by a potential field V(x,y). 
For high-purity materials and at low temperatures, the electron's motion inside the duct is 



a, 



the conductivity of the quantum 



ballistic, and can be described as a scattering process 
circuit being expressible in term of the transmission coefficients of the system. 

Recently, we have developed an S'-matrix approach to stubbed wave guides with de- 
fects, which allows for an accurate numerical solution of the scattering problem even when 
some critical dimension of the system gets large [lOj. We start from the two-dimensional 
Schrodinger equation 

h 2 



2m* 



Vi + V(x,y)}V(x,y) = ES>(x,y) , (0.1) 



where V 2 is the two-dimensional Laplace operator, E the total energy, and m* represents 
the effective mass of the electron in the conduction band. Assuming hard-wall boundary 
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FIG. 1: A stubbed quantum waveguide of width b and infinite length, with a stub of width c and 
length l s . The stub contains a defect with dimensions w X l^. 



conditions, the total wave function ty(x,y) is expanded in terms of the transverse mode 
eigenfunctions in the lead and in the cavity, and the Schrodinger equation (10 .ip is replaced 
by an (in principle) infinite set of coupled, one-dimensional Schrodinger equations. The 
latter can be reduced to linear, algebraic equations matching the wave function and its first 
derivative at the various interfaces delimiting the duct from the cavity, and the empty part 
of the cavity from the region where the potential acts. Thus, the scattering operator for each 
segment in the quantum circuit can be evaluated through linear algebra. The total S'-matrix 
of the device is finally obtained from the partial scattering operators by recursively applying 
the ^-product composition rule, which expresses the overall scattering matrix S in terms of 
the partial scattering matrices and S*^ as 0, lol ] 



where 
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Because of the presence of forward propagating modes only, the evaluation of the scattering 
matrix is numerically stable also for "large" systems. Moreover, the composition rule (j0.2p 
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naturally accommodates a different number of modes in the lead and in the cavity. These 
features are of particular relevance in t he p resent instance, where the stub's width c may 



vary over a rather large range of values 101 ] . 



It is worth to stress here that each block Sy in the scattering operator S is itself a matrix, 
whose elements are labeled by mode or channel indexes. For an incoming wave of unit flux 
impinging from the left, (Sn) nm represents the reflection coefficient towards the left from the 
initial channel m into the final one n, whereas (S2i) nm is the transmission coefficient to the 
right from mode m into mode n. Similarly, (S12) and (S22) nm are the m — > n transmission 
amplitudes to the left and reflection coefficient to the right for an electron incoming from 
the right. Once the transmission coefficients are known, the total conductance G (in units 



2e jh) is given by the Biittiker formula [8! 



9. 



C = £^)l(S 2 iU 2 , (0-4) 

m,n Km 

where and km denote the lead propagation momenta in channel n and m, respectively, 
and the sum is restricted to the open channels in the duct. 

The above S"-matrix approach can be straightforwardly extended to complex energies. 
We used our code to numerically locate the poles of the S"-matrix in the multi-sheeted energy 
surface. In the following, sheets will be specified according to the sign of the imaginary part 

Thus, for a four-channel situation, the 



of the lead momenta in the various channels jl2 |. 
physical sheet, where all the imaginary parts of the momenta are positive, will be denoted 
as (+ + ++), whereas on the sheet ( — h ++) one has Im&f < and Imkf^ > for the other 
three channels. Dominant poles producing resonance effects in the lowest subband, between 
the first and second scattering thresholds E^p and \ are in the fourth quadrant of this 
sheet near the real energy axis, and have < KeE p < E^ . Similarly, dominant poles 

for resonances in the second subband with E^' < E < E^' lie in sheet ( h+) and have 

(2) (3) 

E)p < KeE p < Ej, . We have chosen the value m* = 0.067m e for the effective electron mass, 
which is appropriate for the Al x Gai_ x As/GaAs interface. We verified that convergence is 
attained for both the conductance and the pole positions when four channels are included 
in the external duct, and up to ten channels are taken into account in the cavity. In these 
conditions, the position of the poles in the complex energy plane can be guaranteed with 
an accuracy of the order 10~ 5 . From now on, to exploit the scale invariance of the system, 
all lengths are measured in terms of the waveguide width b, and energies in terms of the 



waveguide fundamental mode = (f ) 2 , and the "tilde" symbol will be used to denote 
adimensional quantities, so that one has for the various thresholds = 1, 2, .... The 
calculations we present refer to a device with l s = 1; a repulsive, double Gaussian defect 

V(x, y) = v oe -P 2 (*-*cf-*Hy-y^ ( 05 ) 

centered in (x c ,y c ) = (0.50,0.25) has been allowed in the cavity. The decay constants along 
the transverse and propagation direction have been fixed at a = 15, (3 = 10, so as to ensure 
that the potential is entirely contained within a region w = 0.3 wide and 1^ = 1 long, 
displaced a distance Y = 0.1 from the lower edge of the guide. The smooth dependence of 
V(x,y) has been taken into account through a slicing technique, i.e., replacing the actual 
interaction with a sequence of pseudodefects having a constant value along the x direction 
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13]. Quite stable results are obtained with N = 10-j-15 slices. In the present calculations 
we have chosen V = 4. 

In Fig. [2] we report the trajectories on the complex energy surface of three S"-matrix 

~ ~ (2) 

poles with varying stub's width c. Pole 1 moves from the upper edge Ej, towards the lower 
edge Ejl of the first subband as c is increased from 1.50 to 5.00. Similarly, pole 2 moves 
downwards from the third threshold passing below the second one as c is increased from 
c = 1.00 up to c = 2.00, whereas pole 3 refer to 1.33 < c < 5.00. In all cases one has the 



"binding" effect typical of an increase of the stub's width [lOj. Note that the three pole 
trajectories appear to be close to each other, but are in fact on different sheets of the energy 
Riemann surface. Pole 1 lies on the ( — h ++) sheet, and can produce resonance effects in 

the first subband, whereas poles 2 and 3 belong to the ( h+) sheet, and are responsible 

of resonance structures in the second subband. As a consequence, pole 2 is a dominant pole 

(2) ~ 

until is passes below Ej, , which happens for c = 1.54; for c > 1.54 it becomes a shadow 

pole, since the ( h+) sheet is more distant from the first subband than the ( — h ++) 

sheet, where the relevant resonance poles may be found. Similarly, pole 3 is shadow for 
c < 2.40, and becomes a dominant pole for greater values of c. In Fig. [2] dominant and 
shadow poles are drawn as full and dashed lines, respectively. 

The change of status of a pole from dominant to shadow pole as it passes a threshold can 
explain the remarkable effects that even small variations of c may have on the conductance 
near threshold. This is illustrated in Fig. [31 where we plot the conductance in the second 
threshold region (3.5 < E < 4.5), in correspondence to c = 1.520, 1.540, 1.541, and 1.560. 



.0.8 I 1 1 1 1 1 1 1 1 1 

1 23456789 10 

Re E 

FIG. 2: Motion of three 5-matrix poles on the Riemann energy surface with varying c. The three 
trajectories correspond to 1.50 < c < 5.00, 1.00 < c < 2.00, and 1.33 < c < 5.00 for pole 1, 2, and 
3, respectively. Note that pole 1 moves on the ( — h ++) sheet, whereas poles 2 and 3 belong to 

the ( h+) sheet. Shadow and dominant poles are drawn as dashed and full lines, respectively. 

The values of c where a pole changes its nature are given in the figure. 

The corresponding conductance profiles are given by the solid, long-dashed, short-dashed, 
and dotted lines, respectively. 

For c = 1.520 pole 2 is dominant, since one has E p ~ 4.14 — 0.392, and produces the 
resonance peak one observes just above threshold. For c = 1.540 and c = 1.541 pole 2 is just 
above (E p ~ 4.002 - 0.276i) and just below (E p ~ 3.990 - 0.270i) threshold. One has that 
the resonance peak is still visible in both cases, which means that the dominant — > shadow 
transition does not prevent the pole from having effects on the observable quantities. For 
c = 1.560 the pole has moved down to E p ~ 3.87 — 0.09i, and it is far away enough from 
the physical region, to have no effects on the conductance, which appears rather flat above 
threshold. It is worth to stress that for these values of c pole 1 is far above the second 
threshold, and cannot influence the conductance profile in the first subband; as a matter of 
fact, in all cases the conductance is practically the same below threshold, and exhibits a cusp 
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FIG. 3: Conductance (in units 2e 2 /h) in the region of the second threshold for c = 1.520 (solid 
line), c = 1.540 (dashed line), c = 1.541 (short-dashed line), and c = 1.560 (dotted line). 

structure, with infinite slope as a function of energy both from above and from below. This 
behavior is indeed discernible at threshold in all calculations, and can be explained much in 
the same way, as one explains threshold phenomena in inelastic scattering processes. When 
a new transverse mode opens up, less energy is available in the propagation direction, so 
that one has the analogue of "endoergic" reactions in inelastic scattering [7j . From Eq. (j0.4p 



one sees that G is linear with respect to the corresponding final momentum kn ■ Since kn 
is related to the total energy E and to the relevant waveguide eigenenergy 6$ by [h]] 

A;« = [2m*(£-a/^ 1/2 , 

one actually expects an infinite derivative of G with respect to E {tJ. 

The effects due to the exchange of role between shadow and dominant poles are illustrated 
in Fig. [31 where we plot G near the second threshold for c = 1.54 (solid line) and c = 1.70 
(dashed line). In the former case one has the resonance peak above threshold due to pole 2, 
as discussed previously; in the latter, pole 2 has moved down to E p ~ 2.68 — 0.022 and has 
no effect on the conductance any longer; pole 1 which moves on the ( — h ++) sheet, on the 
other hand, is now dominant, being located at E p ~ 3.78 — 0.02z, and produces the Fano 



.(0 



dip one observes in Fig. [3j Note that in the first subband one can have the simultaneous 
presence of poles and transmission zeros, which cannot occur when more than a propagating 
mode are active. 
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FIG. 4: Conductance (in units 2e 2 /h) in the region of the second threshold for c = 1.54 (solid line) 
and c = 1.70 (dashed line). 



A similar phenomenon is visible in correspondence to the third threshold. An example 
is given in Fig. [5j where the conductance around the third threshold is plotted for c = 2.37 
and c = 2.42. While a resonance dip is clearly visible in the latter case, no resonance at all 
is discernible for the shorter stub, and only the threshold cusp survives for the conductance 
profile. Such a striking change in correspondence to so small a change in the cavity width 
can be readily explained in terms of a dominant to shadow pole transition. Indeed, for 

c = 2.42 pole 3 of Fig. El is located at E p ~ 8.92 - 0.053i in the ( h+) sheet, and plays 

the role of dominant pole for the second subband. When the stub is shortened, the pole 
moves on its sheet up to E p ~ 9.11 — 0.0232, in correspondence to the third subband, and 
becomes a shadow pole. 

In conclusion, we have demonstrated that the behavior of the conductance near the 
thresholds for the opening of new propagating modes, and its sometimes striking changes in 
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FIG. 5: Conductance (in units 2e 2 /h) in the region of the third threshold for c = 2.37 and c = 2.42. 

correspondence to moderate or even small variations of the stub's width are signals of the 
transition from a dominant to a shadow status of the S'-matrix poles. This result shows that 
concepts and methods of the analytic ^-matrix, widely employed in traditional scattering 
theory, may have their counterpart in the analysis of systems with a confined geometry. 
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